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Abstract. We present results of large scale simulations of the spin-1 Heisenberg 
antiferromagnet on a tetragonal lattice. The stochastic series expansion quantum Monte Carlo 
method is used to calculate equilibrium thermodynamic variables in the presence of an external 
magnetic field. In particular, the low temperature magnetization curve is investigated in the 
quasi-one-dimensional (QID), quasi-two-dimensional (Q2D), and three-dimensional (3D) limits. 
Starting from the 3D limit, the QID (Q2D) limit is achieved by reducing the in-plane (out-of- 
plane) spin coupling strength towards zero. In the QID limit, a Haldane gap appears in the 
magnetization curve at low magnetic field. Additionally, near the saturation field the slope 
of the magnetization curve increases substantially, approaching the infinite-slope behavior of 
a one-dimensional spin-1 chain. A similar (though less pronounced) effect is seen in the Q2D 
limit. We also study the effect of uniaxial single-ion anisotropy on the magnetization curve 
for QID and Q2D systems. Our results will be important in understanding the field-induced 
behavior of a class of low-dimensional Ni-based quantum magnets. 



1. Introduction 

The nature of the ground state phases in interacting quantum many body systems is 
often strongly influenced by the spatial dimensionality. This is a direct consequence of 
enhanced quantum fluctuations that result in novel quantum phases in systems with reduced 
dimensionality that are not present in the isotropic limit. A prominent example of this 
phenomenon is the Haldane gap phase in spin-1 Heisenberg antiferromagnetic (HAFM) chains. 
Haldane proved that the ground state of the one-dimensional spin-1 HAFM is gapped [1] whereas 
there is strong numerical evidence that the ground state in higher dimensions is the gapless Neel 
state [2]. 

In this paper, we study the emergence of the Haldane phase as the quasi-one-dimensional 
(QID) limit is approached, starting from the isotropic spin-1 HAFM and varying the interchain 
coupling continuously. The quasi-two-dimensional (Q2D) limit is also explored by varying the 
inter-layer coupling in a continuous manner. The zero-field ID phase diagram has been well 
mapped |3], and the ID field-dependence has also been reported [1]. The 2D phase diagram is 
also well known |5]. Here, we attempt to extend this knowledge to the QID and Q2D limits, 
and to make a connection to the isotropic 3D limit. In addition to elucidating the nature of the 
emergence of the effects of reduced dimensionality, our results have direct relevance to recent 
experiments on QID and Q2D Ni-based spin-1 Heisenberg quantum antiferromagnets. 
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Figure 1. Magnetization and susceptibility approaching the QID (left) and Q2D (right) limits. 
Results are given for tetragonal lattices of 8 x 8 x 32 (QID) and 16 x 16 x 8 (Q2D). 



2. Theory 

The Hamiltonian for a spin-1 Heisenberg antiferromagnet on an anisotropic cubic lattice is 

where h is the applied magnetic field, Z) is a single-ion anisotropy parameter, and ((u)||) 
represents nearest neighbor spin pairs with in-plane (out-of-plane) bonds. The saturation field 
can be determined by a zero energy cost to flip a single spin in the fully polarized ground state, 
and is given by /i^ = D + 4(2 + Jy). We define a dimensional anisotropy parameter r as 
the ratio J±/J\\ {J\\/Ja.) for QID (Q2D) systems. Thus, in either the QID or Q2D case, the 
low-dimensional limit is achieved as r approaches zero. 

We have performed large-scale numerical calculations of the above Hamiltonian using the 
stochastic series expansion quantum Monte Carlo method with operator-loop updates [6] . In all 
our simulations, either Jy or J_|_ is always set to 1, so that the inverse temperature, /3, is in units 
relative to the isotropic interaction strength (Jy in the QID regime and J±_ in the Q2D regime). 

3. Results 

In Figure [T] we show the magnetization curves as we approach the QID (left panel) and Q2D 
(right panel) limits. Shown in the insets is the uniform static susceptibility, X) which displays 
a peak just before the saturation field. The height of this peak increases as r is diminished, 
the effect being most pronounced for the QID system. Additionally, for the QID system with 
r = 0.01, a Haldane gap appears in the magnetic response at low fields. 

Next, we turn our attention to the low-field behavior of the QID systems, where a spin-gapped 
Haldane phase is expected for ID systems [1]. In Figure [2] we plot the uniform magnetization 
at low field for systems ranging from the isotropic 3D limit to the QID limit. In the extreme 
QID limit (r < 0.01), a spin gap is clearly seen in the magnetic response to applied field. In the 
other limit for r near the isotropic 3D value, no gap is present. Thus, we see a crossover from 




Figure 2. Magnetization at low field for 
the QID system at /3 = 32. A Haldane 
spin-gapped phase emerges for r w 0.02 and 
below. The plateau at m/m^ = 1/16 is due 
to the finite size of our system (4 x 4 x 16). 




Figure 3. Magnetization at low field for the 
QID system with r = 0.01. As the easy-plane 
anisotropy D is increased to 0.40 and above, 
the Haldane gap closes. Results are obtained 
on a 4 X 4 X 16 lattice at /3 = 64. 



3D to ID behavior at a critical value of the anisotropy, Tc, that we estimate to be ~ 0.02, as 
this value separates the two observed limits. Whether or not there is a Haldane gap precisely 
at r = 0.02 is difficult to determine due to possible finite-size and finite-temperature effects. 

Now we determine the effect of adding an easy-plane anisotropy to the zero-field QID system 
in the effective ID limit (r = 0.01). The Haldane gap phase is destroyed as the value of D is 
increased, as shown in Figure [3} Here we note that this occurs rather rapidly, as no gap is seen 
for D > 0.2, when it is known that in ID the gap persists up to D ~ 1 [7]. 

Finally, we show the full field behavior for D > in the QID and Q2D limits in Figure |4j We 
find that the peak in x near the saturation field increases with increasing D for both QID and 
Q2D systems. A broad, shallow peak also emerges in the QID system for D = 2.00. This may 
be associated with the approach of a quantum paramagnetic phase at large D. This spin-gapped 
phase is seen in the QID system at D = 4.00 and in the Q2D system for D = 8.00. 

4. Discussion 

In the Haldane phase at zero field, there is a spin gap A to magnetic excitations. This gap closes 
under an applied field he = A [H]. Below this critical value of the magnetic field the ground 
state magnetization is zero. In ID the Haldane gap is well known to be A = 0.41050(2), in units 
of J|| 13 . Thus he ~ 0.4J|| or, utilizing hs = 4J|| in ID, he ~ O-l/i^. We can estimate he for our 
QID calculations from Figure [2] Taking into account thermal effects, we find he ~ 0.09/is in the 
effective ID limit. This is slightly below, but in general agreement with, the pure ID result. 

We can also compare our estimate of re (~ 0.02, see Figure [2]) with a mean field theory 
(MFT) result of ~ 0.013 [7j. The difference between our value and the MFT calculation can be 
attributed to two likely causes. One is that we are using relatively small system sizes, which 
will tend to show finite size gaps. The other is a tendency for the MFT to overestimate order, 
implying that re from MFT may be a lower bound ^ . 




Figure 4. Magnetization and susceptibility for D > in the QID (left) and Q2D (right) limits. 
Results are obtained on 4 x 4 x 16 (QID) and 8x8x4 (Q2D) lattices at /3 = 16 and r = 0.05. 

5. Conclusions 

We have shown a crossover from 3D to QID physics as the in-plane spin couplings are reduced 
towards zero (relative to the easy-axis spin couplings). The result is the emergence of the 
Haldane phase for r < Tc w 0.02. The behavior of the magnetization curve in the Q2D limit 
was also calculated. The susceptibility peak near the saturation field was found to increase 
with decreasing r. This same effect appears in the QID limit, though much more pronounced. 
We also showed the effect of an easy-plane single-ion anisotropy. In addition to sharpening the 
susceptibility peak, a gapped phase appears for large values of D. 

These results are relevant to understanding the magnetization processes in real materials. In 
particular, the gapless states of the QID and Q2D systems are expected to help characterize the 
magnetic behavior of a class of low-dimensional Ni compounds. We have already seen promising 
qualitative agreement with our Q2D simulations. 
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